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Abstract 

We present an investigation of the propagation of Ultra High Energy Cosmic Rays (UHE- 
CRs) in extragalactic magnetic field (EGMF). We use the IRAS PSCz catalogue in order to 
reconstruct EGMF taking into account power-law dependence between the magnetic field B 
and the infrared luminosity density p^, B = Kpi^. Contrary to the previous works [11] the 
normalization parameters K and (3 have been estimated according to the observable values of 
magnetic field of the galaxy clusters {Bqc = 0.3 — 2 jiG) and assessed values of field in the voids 
of the Local Universe {Bvoid = 10~^^ — 10~^°G). We construct the full-sky maps of expected 
deflection angles of UHE protons with the arrival energies E = 10^° eV and E = 4 ■ 10^^ eV in 
the reconstructed EGMF for two boundary cases of "strong" {Bqc = 2 /iC, Bvoid = 10^^° G) 
and "weak" {Bqc = 0.3 yuG, Bvoid = 10^^^ G) magnetic fields. It is found that average de- 
flections of protons with observable energies above 4 ■ 10^^ eV and maximum energy in sources 
Emax = 10^^ eV are unimportant (comparable with the errors of modern detectors) only in the 
case of "weak" EGMF model. 



1 Introduction 

The nature of UHECRs with E > 10^^ eV is one of the most challenging problems of modern 
astrophysics [1, 2]. At these energies the cosmic ray particles strongly interact with the cosmic 
microwave background and their flux should be severely attenuated according to GZK cutoff 
of UHECR energy spectrum [3, 4]. AGASA detector in Japan has registered 11 events with 
energies above 10^*^ eV [5]. If galactic and extragalactic magnetic fields are not enough strong, 
then protons of these energies may not be significantly defiected so that some of the cosmic rays 
may point back to their sources. Therefore, the EGMF is of great importance for the search of 
UHECRs sources. 

A few different methods of EGMF reconstruction are presented in literature. The authors 
[6, 7] described the magnetic field distribution by simple models that were connected with such 
large-scale structures as the Local Supercluster, the Supergalactic wall etc. Recently, several 
groups started developing physically more realistic models based on numerical simulations, 
combining the magneto-hydrodynamics of the intergalactic medium with N-body simulations 
of the driving gravitational dynamics of the dark matter [8, 9, 10]. The method that was 
presented in work [11] is realistic and relatively simple simultaneously. In this work the EGMF 
simulation was carried out on the basis of power-law dependence between the galaxy infrared 
luminosity density and strength of the magnetic field. 

We use this approach in our work to reproduce the EGMF, but we realize a calibration of 
the field energy other way. The authors [11] set a normalization of magnetic field B = 0.4 /xG 
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in a 1 Mpc^ cube in the centre of the Virgo cluster and took the parameter j3 = 2/3 in the 
dependence between the infrared luminosity density and EGMF . We cahbrated the simulated 
EGMF using observable magnetic field values in the centre of the Coma galaxy cluster and 
assessed values of the field in the voids of the Local Universe. 



The EGMF are little known theoretically and is difficult for observations. The results of [9, 12] 
show that the magnetic field traces baryon density. On the other hand, the Faraday rotation 
measurements of polarized radio sources placed within cluster of galaxies and observations of 
hard X-ray emission from clusters of galaxies points at connection between magnetic field and 
distribution of galaxies [13]. We have reconstructed the EGMF, from the hypothesis about its 
origin as a result of magnetized plasma injection in the cxtragalactic medium (within 1 Mpc 
aroung the galaxy) via supernova explosions and galactic wind action [14]. If the galactic 
winds from supernova explosions magnetize the intergalactic medium in the same way as they 
magnetize the interstellar medium, then, under the assumption of equipartition, the magnetic 
field in the intergalactic medium should have the magnetic field energy density comparable 
to the energy density of cosmic rays [15, 16]. The authors [11] assume that the luminosity 
density is proportional to the gas density. In our work we used the IRAS PSCz catalogue 
within 100 Mpc [18] for construction of EGMF model similar to [11]. We applied power-law 
dependence between magnetic field B and infrared luminosity density pl, at a point "r*, in the 
form of: 



where plo and Bq is the values of infrared luminosity density and magnetic field, and P is the 
free parameter. For calibration of this equation we took the conditions in the Coma cluster 
centre (/ = 58°, b = 88°, r = 96 Mpc, plo = 1.1 • Mpc^) and in the voids of the 

Local Universe. By comparing the Coma radio synchrotron spectrum with excess radiation, 
interpreted as the inverse Compton scattering of cosmic microwave background photons by 
same electron population, volume-averaged magnetic field of i? ~ 0.1 pG have been deduced 
in [20], 0.2 pG in [21], and O.ApG in [22]. Faraday rotation measurements gave B r^2pG [17], 
8.5 pG [19]. The obtained values of magnetic field may be divided into two groups: B < IpG 
and B > IpG. We took the value Bq — 0.3 pG as representative for the first group, and 
Bq = 2.0 pG for the second group. 

The data about Local Universe voids were taken from [25]. The mean value of void infrared 
luminosity density pvoid is equal to 9.9 • IO^Lq/Mpc^. We took the value Bvoid = 10^^° G 
as representative for the upper hmit of EGMF in voids and the lower limit Bvoid — 10~^^ G 
for the magnetic fields in the voids from [24]. Taking into account the data described above 
the two approximate dependences between magnetic field and infrared luminosity density were 
obtained: the model of "weak" field {Bq = 0.3 pG, Bvoid = lO'^^G) - 



2 Models of the extragalactic magnetic field 
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Figure 1: The magnetic field distribution along the three fiducial lines through the Virgo cluster, 
the Perseus cluster and the Coma cluster from [11] (left panel) and from our simulation (right 
panel). The shaded and unshaded regions on the right panel represent the models of the "weak" 
and "strong" magnetic fields, respectively. 



and the model of "strong" field {Bq = 2.0 /xG, Bvoid = 10"^° G) - 

B{r) = 2.o ( '4^^ — X\g, 

The magnetic field distribution along the three fiducial fines through the Virgo cluster, the 
Perseus cluster and the Coma cluster are shown in the Fig. 1. To avoid zero magnetic field 
value we assigned the minimal EGMF {Byoid) to volumes with Pl(^) = 0. 

It is reasonable that direction of the EGMF induction is random. Coherent lengths of 
EGMF Ic in the galaxy cluster regions lie in the range from 10 to lOO/cpc [26, 27], and are 
approximately equal to 1 Mpc for intergalactic medium with B < 10~^G [26]. We have used 
the following conditions: if _B < 10~^ G then value of Ic is equal to 1 Mpc, \l B > 10"^ G then 
Ic took on the values 250, 50 and lOkpc. 



3 The deflection of UHECRs in the extragalactic mag- 
netic field 

We have studied the UHECR deflections in the cases of different EGMF models. For this we 
injected antiprotons from the Earth and followed their trajectory taking into account energy 
losses of the UHECRs due to the interactions with cosmic microwave background photons [23] . 
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Figure 2: The deflection angles of UHE protons with arrival energy E — 4-10^^ eV (top panel) 
and E — 10^° eV (bottom panel) through propagation in "weak" model of extragalactic mag- 
netic field with Ic — 50 kpc. Both maps are smoothed with Gaussian kernel of 3°. The coordinate 
system is galactic. Deflections in degrees. 

Table 1: The average deflection of the UHECRs in the different EGMF models 







Energy, cV 


Ic, kpc 


Model 


4 • 10^9 


6 • W 




3 ■ 10^0 


250 


weak 


9.5° 


7.1° 


2.5° 


0.4° 


strong 


23° 


19° 


11° 


2.9° 


50 


weak 


5.3° 


3.4° 


1.2° 


0.2° 


strong 


18° 


14° 


7° 


1.3° 


10 


weak 


2.4° 


1.5° 


0.5° 


0.1° 


strong 


12° 


9.2° 


3.3° 


0.6° 
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The trajectories were followed until the distances of particles from the Galaxy reached 
100 Mpc or their energies reached 10^^ eV. The average deflections of the UHECRs significantly 
depend on the coherent length value Ic and type of EGMF model (Tab. 1). So, the average 
deflection for the "strong" model is 3 - 6 times as much as that for the "weak" model. The 
deflection maps of protons in the case of the "weak" magnetic field model for l^. = 50kpc are 
presented in the Fig. 2. As can be seen from these maps, regions with significant deflections 
(above 15° for UHECRs with arrival energy 4 ■ 10^^ eV and above 7° for 10^° eV) correspond to 
locations of main galaxy clusters and near galaxies with high infrared luminosities. Although 
the most part of the sky area is occupied by regions with deflection values less then 4° for 
E = A - 10^^ eV and 1° for E = 10^° eV. Thus it is appropriate to search of correlation between 
the UHECRs (with arrival energy above 4 • 10^^ eV) arrival directions and their possible sources 
in the case of "weak" model. 

4 Conclusion 

In this work we have investigated the influence of EGMF on the UHECRs propagation. For 
this purpose we have reconstructed EGMF taking into account power-law dependence between 
magnetic field and infrared luminosity density of galaxies in the Local Universe. As distinct 
from [11] this dependence was normalized to the observable values of the magnetic field in the 
galaxy clusters and the assessed values of the field in the voids of the Local Universe. We have 
found the value of exponent /3 in power-low dependence for the two boundary cases of magnetic 
fields in the Coma cluster and in the voids. For the "weak" EGMF model {Bgc = 0.3/xG, 
Bvoid = 10~^^ G) the parameter f3 is equal to 1.8 and for the "strong" field model {Bgc = 
2.0/iG, Bvoid = 10~^° G) P = 1.4. The average values of UHECR deflection depend on the 
type of EGMF model. It was shown, that deflection angles for protons with the arrival energy 
equal to 4 • 10^^ eV signiflcantly exceed errors of modern detectors in case of "strong" EGMF 
model. For "weak" EGMF model the deflections lie mainly in the range of 1° — 5°, therefore 
charged particle astronomy should be possible in that case. 

On the basis of the average deflection comparison we can conclude, that our "weak" model 
of EGMF agrees with the model that is presented in [11]. On the other hand our results differ 
from [9], where the authors constructed EGMF model using simulations of large-scale structure 
formation to study the build-up of magnetic flelds in the intergalactic medium. In average our 
deflection angles are larger then obtained in [9]. Our results are in good agreement with [10, 
12] where the signiflcant UHECR deflections in simulated EGMF were found. 
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